Neutron capture cross section of unstable ^^Ni: implications for stelleir nucleosynthesis 
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The ®^Ni(n, 7) cross section has been measured for the first time at the neutron time-of-flight 
facihty n_TOF at CERN from thermal neutron energies up to 200 keV. In total, capture kernels 
of 12 (new) resonances were determined. Maxwellian Averaged Cross Sections were calculated 

for thermal energies from kT = 5 kcV to 100 kcV with uncertainties around 20%. Stellar model 
calculations for a 25 star show that the new data have a significant effect on the s-process 
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production of ^^Cu, ^"Ni, and '^^TjU in massive stars, allowing stronger constraints on the Cu yields 
from explosive nucleosynthesis in the subsequent supernova. 

PACS numbers: 25.40.Lw, 25.40.Ny, 26.20.Kn, 27.50.+C, 28.20.Fc, 29.30.Hs 



The weak component of the astrophysical s process 
observed in the Solar System abundance distribution in- 
cludes the s-process species between Fe and Sr (60 < 
A < 90) |I] . Most of them are generated in massive stars 
during convective core He-burning and convective shell 
C-burning via the activation of the neutron source re- 
action ^^Ne(a,n)^^Mg The long-lived radioisotope 
isotope "^^Ni (^1/2=101.2 ± 1.5 yr [7]) is located along 
the neutron capture path, and in typical weak s-process 
conditions it may become a branching point, when the 
neutron capture timescale is comparable with its stellar 
/3-decay rate. 

In particular, at the end of He core burning the neu- 
tron source ^^Ne(a,n)^^Mg is activated at temperatures 
around 0.3 GK (GK = 10^ K), corresponding to a 
Maxwellian neutron energy distribution for a thermal 
energy of kT = 26 keV. At this stage, neutron densi- 
ties are too weak for a subsequent neutron capture on 
^''Ni (with central peak neutron density in the order of 
10'' cm-3, e.g. [g [6]) and more than 90% of the *^^Ni 
produced decays to ^^Cu. However, the s-process mate- 
rial is partly reprocessed during C shell burning, where 
the ^^Ne(a, n)^^Mg neutron source is reactivated at much 
higher temperatures of about 1 GK, corresponding to a 
thermal energy of kT = 90 keV. 

During this second stage, neutron densities are orders of 
magnitudes higher, reaching a maximum of 10^^"^^ cm^"^ 
[8,. At the ^^Ni branching point, the high neutron den- 
sities favor neutron capture producing ^"^Ni, bypassing 
the production of ^'^Cu despite of the strong temperature 
dependence of the ^'^Ni /3-decay rate (at C shell temper- 
atures the half-life of ^'^Ni decreases to few years [H]). 
In these conditions, the amount of ^^Cu generated dur- 
ing He core burning is partially depleted in the C shell, 
but the final ^'^Cu abundance will increase thanks to the 
later radiogenic decay of the ^'^Ni accumulated in the C 
shell burning phase [8]. In Fig. 1, we show the neutron 
capture path in the Ni-Cu-Zn region during He core and 
at high neutron density during C shell burning. Up to 
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FIG. 1. The s-process reaction path in the Ni-Cu-Zn region 
during He core burning (dashed lines) and C shell burning 
(solid lines). 



now the stellar cross section of ^■^Ni(n, 7)^''Ni relied on 
calculations or extrapolations of experimental values at 
thermal neutron energies (0.025 eV) [TDHH]- Theoretical 
predictions for the Maxwellian Averaged Cross Section 
(MACS) at kT = 30 keV are ranging from 24 to 54 mb 
|13H17j . The currently recommended value quoted by the 
compilation KADoNiS [TS] is 31 ± 6 mb. Because such 
calculations are vulnerable to large systematic uncertain- 
ties, measurements have been attempted at Los Alamos 
National Laboratory [TH] and at CERN. In this letter we 
report on the first experimental results for the ®^Ni cross 
section at stellar energies obtained at the n_TOF facility 
at CERN. 

The measurement was performed at the neutron time of 
flight facility n_TOF located at CERN. Neutrons are pro- 
duced via spallation reactions of 20 GeV/c protons from 
the Proton Synchrotron with a massive Pb target. With 
the high intensity of the pulsed proton beam, the repe- 
tition rate of 0.4 Hz, a short proton pulse width of 6 ns, 
and a neutron flight path of 185 m, the n_TOF facility is 
unique for the combination of a high neutron energy res- 
olution and a high instantaneous neutron flux. A more 
detailed description of the facility can be found in Ref. 
[20] and references therein. 

The ^'^Ni sample was produced by irradiating highly en- 
riched ^^Ni in a thermal reactor [HI [HI [22]. Since the 
irradiation took place more than 20 years ago, the orig- 
inal ^'^Ni fraction had partially decayed to ^^Cu. To 
avoid background due to ^^Cu(n, 7)^^Cu reactions, this 
^^Cu impurity has been chemically separated prior to 
the (n, 7) measurement. The originally metallic target 
was dissolved in concentrated nitric acid and the copper 
fraction was precipitated as CuS using gaseous H2S. The 
remaining solution was treated with NaOH to precipi- 
tate Ni(0H)2, which was calcinated at 800 °C to form 
NiO. By means of mass spectrometry, the ^^Ni/^^Ni ra- 
tio in the sample was determined to 0.123 ± 0.001%, and 
the contribution from other Ni isotopes was found to be 
< 1%. In total, 1156 mg NiO powder were encapsulated 
in a thin-walled cylinder made of PEEK (Polyetherether- 
ketone, net mass 180 mg) to produce a sample 20 mm in 
diameter and 2.2 mm in thickness. 

The neutron capture yield was measured as a function of 
neutron energy by detecting the prompt capture 7 rays 
with a pair of liquid CgDg scintillation detectors. These 
detectors are optimized to exhibit a very low sensitivity 
to neutrons, thus minimizing the background produced 
by neutrons scattered on the sample [23]. The depen- 
dence of the detection efficiency on 7-ray energy and the 
effect of the 7-ray threshold of 250 keV were corrected 
using the Pulse Height Weighting technique [211 [25] ■ By 
application of a pulse-height dependent weight on the de- 
posited 7 energy the detection efficiency becomes a linear 
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function of the excitation energy of the compound nu- 
cleus, £ K, k X Ec- Choosing k = 1 MeV^^, the capture 
yield can be obtained as 



where Cw are the weighted, background-subtracted 
counts, $ denotes the relative neutron flux, and N a 
normalization factor for the absolute capture yield. The 
normalization factor was determined via the saturated 
resonance technique |26| in an additional run with a Au 
sample of the same size as the Ni sample. The Au sample 
was chosen such that the gold resonance at 4.9 eV is sat- 
urated, which means that all neutrons of that energy are 
absorbed in the sample, thus providing a measure for the 
absolute neutron flux at 4.9 eV. The energy dependence 
of the neutron flux was measured relative to the stan- 
dard cross sections ^'^B(n,a) and ^Li(n, a) up to 150 keV 
and ^^^U(n, /) at higher energies. Because the size of 
the neutron beam widens slightly with neutron energy, 
the normalization factor N, related to the fraction of the 
neutron beam intercepted by the sample, changes as well. 
This effect was taken into account by simulations of the 
neutron beam profile [20 . 

The experimental background was determined in dedi- 
cated runs with an empty PEEK container, with a ^^Ni 
sample of the same diameter, and in runs without neu- 
tron beam. Additionally, the neutron capture yield has 
been measured with a set of neutron filters located about 
50 m upstream of the sample. These W, Mo, and Al fil- 
ters are thick enough to exhibit black resonances at cer- 
tain energies, so that all neutrons in these windows are 
completely removed from the beam and do not reach the 
sample. Accordingly, the level in the corresponding dips 
in these spectra is expected to represent the experimen- 
tal background. When the run with filters was repeated 
using only the empty container, the same background 
level was observed in the filter dips, thus confirming the 
background measured with the empty sample container. 
Figure [2] shows a comparison of the capture yield of the 
^■^Ni sample, the empty PEEK container and the ^^Ni 
sample. The background from the radioactivity of the 
^^Ni sample and from ambient radiation, which was ob- 
tained from a measurement without neutron beam, is 
given as well. Between 100 eV and 2 keV four resonances 
are visible in the spectrum of the ^'^Ni sample, which 
are obviously not correlated with the ^^Ni(n, 7) content. 
The first of these resonances (marked by an arrow) can 
be attributed to a resonance in the ^^Ni(n, 7) reaction 
|27j . expected at that neutron energy and compatible 
with the measured 0.03% impurity of ^^Ni in our sample. 
The other three resonances are clearly attributable to the 
^'^Ni(n, 7) channel. This holds also for several other res- 
onances up to neutron energies of 55 keV, for which the 
capture kernels 
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FIG. 2. (Color online) Capture yield of the ''^Ni sample 
(black) compared to the empty sample holder and to the spec- 
trum obtained with a pure ''^Ni sample. The background from 
the activity of the sample and from ambient radiation (blue) 
is almost negligible at keV energies. The capture yield of the 
^^Ni sample was scaled for the ®^Ni mass present in the ^'^Ni 
sample. The first resonance at 203 eV (marked by an arrow) 
is assigned to a small impurity of ^^Ni in the sample. 

characterizing the strength of the resonance, could be 
deduced by a resonance shape analysis (RSA) with the 
R-matrix code SAMMY [28j (Table |l]). The capture ker- 
nel is determined by the spin statistical factor t/g, the 
neutron width r„, and the radiative width F..^. For two 
resonances also the orbital angular momentum derived 
from the shape of the resonance, is given. The neutron 
energy interval between 2 and 8 keV is dominated by 
the strong resonance in ^^Ni(n, 7) at 4.6 keV, therefore 
smaller resonances in ^^Ni(n,7) might be invisible due 
to this background. In summary, 12 levels in ^^Ki were 
observed for the first time. 

TABLE I. Resonance energies Er (laboratory energy) and 
capture kernels for the ®''Ni(n, 7) reaction. For resonances 
marked with an asterisk the orbital angular momemtum £ = 
could be deduced from the resonance shape. 



Er (CV) 


(meV) 


Er 


A^ (meV) 


397.96 ± 0.04 


5.7 ±0.4 


9776 ± 3 


100 ± 10 


587.25 ±0.09* 


340 ± 20 


13984 ± 3 


131 ± 45 


1366 ± 1* 


810 ± 40 


17127 ±4 


108 ± 59 


8634 ± 2 


45 ±9 


19561 ±6 


130 ± 20 


8981 ± 3 


50± 10 


32330 ± 10 


500 ± 200 


9154 ±4 


43 ±9 


54750 ± 30 


700 ± 200 



As a consequence of the small sample mass, the signal 
to background ratio starts to deteriorate already above 
10 keV. Accordingly, it is increasingly difficult to iden- 
tify resonances with confidence at higher energies. Thus, 
MACSs were calculated using resonance parameters only 
below 10 keV, whereas averaged cross section data have 
been determined from 10 keV to 200 keV. These data 
were obtained by subtraction of the yield measured with 
the ^•^Ni sample after it had been properly scaled for the 
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^^Ni content of the ^"^Ni sample. The background due to 
oxygen is neghgibly small because of its very small (n, 7) 
cross section. 



TABLE II. Maxwellian Averaged Cross Sections (in mb) of 
®'^Ni(n, 7) compared to previously recommended values based 
on theoretical predictions [TS]. The respective contributions 
from resonances below 10 keV (RC) are listed separately. Un- 
certainties are 1 a. 



kT 
(keV) 


KADoNiS 




This Work 


RC 


Total 


5 


112 


174 ±6 


224 ± 8stat ± 45syB 


10 


66 


51 ± 2 


129.5 ±7.1 ±25.9 


15 


50 


24 ± 1 


101.3 ±6.9 ±20.3 


20 


41 


14 ± 1 


85.5 ± 6.4 ± 17.1 


25 


35 


9.3 ±0.4 


74.9 ±5.9 ±15.0 


30 


31±6 


6.6 ±0.3 


66.7 ±5.4 ±13.3 


40 


25 


3.8 ±0.2 


54.5 ± 4.6 ± 10.9 


50 


20 


2.4 ±0.1 


45.6 ±3.9 ±9.1 


60 


17 


1.7 ±0.1 


38.8 ±3.4 ±7.8 


80 


13 


0.97 ±0.05 


29.1 ± 2.7 ± 5.8 


100 


10 


0.63 ±0.03 


22.5 ±2.1 ±4.5 



The MACSs for thermal energies from kT = 5 to 
90 keV are listed in Table [ll] together with the theoret- 
ical predictions in the KADoNiS compilation [18]. Our 
results are approximately a factor of 2 higher than the 
calculated cross section. The total systematic uncertain- 
ties in our results of 20% are mainly due to subtraction 
of the background and the effect of sample impurities - 
particularly in the region between 2 and 8 keV, where 
the spectrum is dominated by the 4.6 keV resonance in 
^^Ni(n,7). Comparably minor contributions to the sys- 
tematic uncertainty are caused by the neutron flux (3- 
5%), the Pulse Height Weighting technique (2%), the flux 
normalization (1%), and the ^^Ni/^SNi ratio (1.6%). 

The impact of our new results was investigated for 
the s process in a full stellar model for a 25 Mq star 
with an initial metal content Z =0.02. The complete 
nucleosynthesis was followed with the post-processing 
NuGrid code MPPNP HHj- The steUar rates were 
obtained by combining the measured ^•^ssjsjj^^^^ ^^64]\jj 
cross sections and theoretically predicted contributions 
to the stellar rate due to ^■^Ni*(n, 7)^**^ reactions 
as described in Ref. While the contribution of 

63gs]\jjj-j^^ ^■j64]v,jj reactions to the stellar rate is still 
around 90% at He Core burning temperatures, it drops 
to around 40% at the higher temperature in the C 
Shell burning phase. Because of the larger uncertainties 
of the ^■^Ni*(n, 7)^'*Ni cross sections, the uncertainty 
of the stellar rate increases with temperature. Apart 
from the reaction cross sections, the final abundance 
pattern is also affected by the temperature dependence 
of the radioactive decay rates under stellar conditions. 
In the investigated mass region the concerned rates for 
the /3--decay of ^^Ni and the /3+//3--decays of '^''Cu 
have been adopted from Ref. [9]. By variation within 
reasonable limits [8] it was found that the decay rates 
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FIG. 3. (Color online) (a) Final isotopic s-process distribu- 
tions using the new measured '''^Ni MACS (red circles) and 
the MACS quoted by KADoNiS (blue squares) [18]. The dis- 
tribution is normalized to solar system abundances, (b) Ratio 
of the two distributions in Panel (a) , zoomed in the Ni-Cu-Zn 
mass region. Isotopes of the same element are connected by 
solid lines. 



of both isotopes have a comparably small effect on the 
investigated abundances, because the reaction flow in 
the ®^Ni branching is governed by the neutron density 
conditions, which lead either to much lower or much 
higher (n, 7) rates during core He and shell C burning, 
respectively. 

The calculated abundance distribution from Fe to Zr 
shown in Fig. [3| represents the s abundances after core 
He and shell C burning, i.e. prior to the supernova 
explosion, at a point where the nucleosynthesis yields are 
well characterized by the model [5]. The distribution is 
compared in Fig. [3] to the one obtained with the neutron 
capture rates of the KADoNiS evaluation [18] . The ratio 
of the two distributions in the lower panel of Fig. [3] 
shows that the new ^^^Ni cross section affects only a few 
isotopes between Ni and Zn. An enhancement of about 
20% is found for ^'^Ni, while ^'^Cu is depleted by about 
15%. As the ^^Cu yields remain essentially unchanged, 
the isotopic ratio ^^Cu/^^Cu is correspondingly reduced 
at the end of shell C burning. ^'*Zn is depleted as well 
(by about 30%), because ^'^Cu and ^^Zn are populated 
by the nucleosynthesis channel following the /3~ branch 
62Ni(n,7)63p^i(/3-)63Cu(n,7)64Cu(/3-)^'*Zn. However, 
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the s-process contribution to Zn remains marginal as 
this isotope results predominantly from later explosive 
nucleosynthesis during Core Collapse Supernovae [5T] . 
Also the propagation effect of the new MACS of ^'^Ni on 
heavier s-process species is rather small, of the order of 
a few percent. 

Although the s process component at the end of 
convective shell C burning is well defined by these cal- 
culations, the abundances in the Ni-Cu-Zn region may 
be affected by following burning stages (for instance the 
possible merging of shells [35]) and by the subsequent 
supernova explosion before enriching the interstellar 
medium. Given the complexity of this scenario, the final 
abundances are yet subject to considerable uncertainty 
as emphasized by several sensitivity studies [SI 133] ■ 
Nevertheless, the present results represent a fundamental 
improvement in constraining the weak s-process compo- 
nent from the convective core He burning and convective 
C shell burning phases. A better knowledge of the 
pre-explosive weak s-process component will allow to 
better define also the following explosive contribution to 
the copper inventory, once robust theoretical predictions 
are compared with spectroscopic observations. Another 
relevant observational constraint is given by the copper 
isotopic ratio in the Solar System, where s process in 
massive stars provide the dominant contribution 
and references therein). 



In summary, we measured the energy-dependent 
^^Ni(n,7) cross section at the n_TOF facility providing 
the first experimental results for MACSs at stellar 
neutron energies. The MACSs ranging from kT = 5 to 
90 keV exhibit total uncertainties of 20 — 22% and are 
about a factor of 2 higher than the theoretical prediction 
of the KADoNiS compilation. Our results improve one 
of the main nuclear uncertainties affecting theoretical 
predictions for the abundances of ^'^Cu, ^^Ni and ^^Zn 
in s-process rich ejecta of core collapse supernovae. 
Furthermore, these results are a fundamental step to 
constrain the contribution from explosive nucleosynthe- 
sis to these species. 
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